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Abstract 

The recent analysis of nuclear distortions in DIS off nuclei revealed a breaking of the 
conventional hard factorization for multijet observable. The related pQCD analysis 
of distortion effects for jet production in nucleus-nucleus collisions is as yet lacking. 
As a testing ground for such an analysis we consider the Abelian problem of higher 
order Coulomb distortions of the spectrum of lepton pairs produced in peripheral 
nuclear collisions. We report an explicit calculation of the contribution to the lepton 
pair production in the collision of two photons from one nucleus with two photons 
from the other nucleus, 27 + 27 —* 1 + 1~. The dependence of this amplitude on the 
transverse momenta has a highly nontrivial form the origin of which can be traced to 
the mismatch of the conservation of the Sudakov components for the momentum of 
leptons in the Coulomb field of the oppositely moving nuclei. The result suggests that 
the familiar eikonalization of Coulomb distortions breaks down for the oppositely 
moving Coulomb centers, which is bad news from the point of view of extensions 
to the pQCD treatment of jet production in nuclear collisions. On the other hand, 
we notice that the amplitude for the 27 + 27 — > 1 + 1~ process has a logarithmic 
enhancement for the lepton pairs with large transverse momentum, which is absent 
for 717 + m l ~ ¥ processes with m,n > 2. 

We discuss the general structure of multiple exchanges and show how to deal with 
higher order terms which cannot be eikonalized. 
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I. INTRODUCTION 

The exact theory of Coulomb distortions of the spectrum of ultrarelativistic lepton pairs 
photoproduced in the Coulomb field of the nucleus has been developed by Bethe and Max- 
imon It is based on the description of leptons by exact solutions of the Dirac equation 
in the Coulomb field (see e.g., the textbook j2|). In the Feynman diagram language one has 
to sum multiphoton exchanges between produced electrons and positrons and the target nu- 
cleus. For ultrarelativistic leptons this reduces to the eikonal factors in the impact parameter 
representation. In the momentum space the same eikonal form leads to simple recurrence 
relations between the (n + 1) and n-photon exchange amplitudes [2|, the incoming photon 
can be either real or virtual. There are two fundamental points behind these simple results: 

i) The lightcone momenta of ultrarelativistic leptons are conserved in multiple scattering 
process (i.e., if the nucleus moves along the n_ lightcone and the produced leptons 
move along the n + lightcone, then the p + components of the lepton momenta are 
conserved) . 

ii) The s-channel helicity of leptons is conserved in high energy QED (see the textbook 
P). It is the last property by which distortions reduce to a simple eikonal factor. 

The same properties allow one to cast the pair production cross section in the dipole 
representation [4]. They have also been behind the color dipole pQCD analysis of nuclear 
distortions and the derivation of nonlinear /c^-factorization for multijet hard processes in 
DIS off nuclei J 

As was shown in jsj, in certain cases of practical interest the so-called abelianization 
takes place. Specifically, the hard dijet production in hadron-nucleus collision is dominated 
by a hard collision of an isolated parton from the beam hadron simultaneously with many 
gluons from the nucleus which belong to different nucleons of a target nucleus. None the 
less, at least for the single-particle spectra, the interaction with a large number of nuclear 
gluons can be reduced to that with a single gluon from the collective gluon field of a nu- 
cleus, i. e., the nonlinear k i -factorization reduces to the linear one and in terms of the 
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collective glue one only needs to evaluate the familiar Born cross sections. The extension 
of nonlinear /c_|_ -factorization for hard processes from hadron-nucleus collisions to collisions 
of ultrarelativistic nuclei is a formidable task which has not been properly addressed so far. 
The lightcone QED and QCD share many properties, and here we address a much simpler, 
Abelian, problem of Coulomb distortions of lepton pairs produced in peripheral collisions of 
relativistic nuclei. 

The process of lepton pair production in the Coulomb fields of two colliding ultrarelativis- 
tic heavy ions was intensively investigated in the past years Such 
activity is mainly connected with new possibilities opened with operation of such facilities 
as RHIC and LHC. Despite the high activity in this area the issue of correct allowance for 
final state interaction of produced leptons with the colliding ion Coulomb fields is lacking 
yet. The main results obtained so far in this direction are the following: 

i) The produced lepton pair interacts with the Coulomb fields of the ion and the cor- 
responding corrections have a noticeable impact on the cross section of the process 
under consideration at finite energies [lo| . 

ii) The perturbation series corresponding to multiple interaction of a produced pair with 
Coulomb fields can be summed and the result can be cast in the eikonal-like form [ 14] , 

n 

if one restricts ourself to terms growing with energy in the cross section [l^. In QED 
such an approximation can be considered as satisfactory one, but it does not work in 
QCD and the problem of higher order corrections in pair production demands further 
investigation. 



In our paper [121] . we cited the amplitude M$ which is irrelevant in leading and next-to 
leading logarithmic approximations in QED. Nevertheless, the knowledge of such a kind 
contributions becomes important for similar processes in QCD with multigluon exchanges 
between the color constituents of each of the colliding hadrons and the created quark- 
antiquark pair. Thus, the main motivation of the present paper is a further investigation of 
multiple exchanges and their impact on the lepton pair yield in the ultrarelativistic heavy ion 
collisions, an issue which is useful not only in understanding the electromagnetic processes, 
but has a wide application in QCD. 

We did not consider the case when one of the ions radiates a single photon and other 
one radiates an arbitrary number of photons absorbed by a created pair The photon 
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exchanges between the ions also were not taken into account 

Our paper is organized as follows. In SecjnJ we consider the case when each of the 
colliding ions radiated two photons which created the lepton pair. We derived the relevant 

(2) 

amplitude M)J using the powerful Sudakov technique well suited for calculations of the 
processes at high energies. 

In Sec JIH| we studied the wide-angle limit in pair production kinematics corresponding 
to the case of large transverse momenta of pair components. In these limits the results 
are much more transparent than in the general be seen from the form of the 

differential cross section which is also presented. 

In SeciIV| we discuss the generalization of the process under consideration to the case, 
when the number of exchanged photons by each ion exceeds two. 



II. THE LEPTON PAIR PRODUCTION 

We are interested in the process of lepton pair production in the collision of two relativistic 
nuclei A, B with charge numbers Z\, Z2 

A( Pl ) + B(p 2 ) -> r(g_) + I+(g + ) + A(p[) + B(p' 2 ), (1) 

with kinematical invariants 

s = (pi+ P2) 2 , ql = (Pi - Pif, ql = (P2 - p' 2 ) 2 
Sl = (q + + q-)\ p 2 1= p? = Ml p\= v 'l = Ml, ql = m\ (2) 

We are interested in peripheral kinematics, i. e., 

s Ml Ml \q\l \qll^> m 2 (3) 

which corresponds to small scattering angles of ions A and B. 

It is convenient to use the Sudakov parameterization for all 4-momenta entering the 
process (JTJ) 

qi = aip2 + bipi + gix, <?2 = a2p2 + b 2 pi + <?2±, 
ki = a x p 2 + Pipi + ku_, k 2 = a 2 p2 + P2P1 + k 2 ±, 

q± = a±p 2 + P±pi + q±±, (4) 
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with lightcone 4-vectors obeying the conditions 



v\=v\ = 0, Pi,2.g± = 0, 2pi.p2 = s. 



A. The pair production by 4 photons 

Let us consider the creation of the lepton pair by four virtual photons (Fig. [TJ. The 
photons with momenta ki, q\ — k\ (in the latter article, referred to as photons 1 and 2) are 
emitted by the ion A and the photons with momenta ki-, q 2 — k 2 (referred as the photons 3 
and 4) by the ion B. The main contribution to the cross section gives the following regions 
of the Sudakov variables: 

ai</3i~6i, /?+ + /3_ = &i, 

(3 2 < «2 ~ «2, «+ + «_ = a 2 , (5) 

|a 1 |<a 2 , |6a|<6i, Qi± = Qi, <?i + <? 2 = 9+ + 9-, 

Q4- 9 2 
a± = —5-, <? ± > m . 

sf3± 

Hereinafter q i denotes the 2-dimensional transverse part of any considered momenta. For 
defmitness, we suggest @ + , /3_ > 0, which corresponds to the situation when the pair moves 
along the ion A (the momentum p±). Bearing in mind a possible extension to pQCD we 
neglect the lepton masses whenever appropriate. 

The contribution to the matrix element of such a set of Feynman diagrams (FD) reads 

(2) . (z 1 z 2 ) 2 (4 TO ) 4 r d%d% 

(2) 15 (2vr)8 J kfkl( qi - hy(q 2 - k 2 y 

x -u^O^ip^ip'^ (6) 

To see the proportionality of the matrix element © to invariant energy s, we use the Gribov 
representation for virtual photon Green functions 

/2\ 4 

# W iflWi5Wi#™i ~ ("J Pl l xPluPlp l Pla 1 P2,xxP2v 1 P2pP2c- (J) 

Numerators of Green functions of the nuclei A can be written as s 2 Ni with N\ = 

u v (p' l )p 2 u ri (p\) / s , l-^i | 2 = 2 and a similar expression takes place for the nuclei B. The 

v 
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denominators of virtual photon Green functions in the considered kinematics depend only 
on transverse components of the corresponding 4-vectors, thus 

k\k\{ qi - k l )\q 2 - k 2 f = k\kl{ qi - k x f{q 2 - k 2 f. 

(2) 

There are 24 FD contributing to M}J. Instead of them it is convenient to consider 24*2*2 = 
96 FD which take as well the permutations of emission and absorption points of exchanged 
photons to the nuclei (Fig. EJ). Then the result must be divided by (2!) 2 . This trick 
provides the convergence of integrals over f3 2 



oo 

1 



5S /dA 



+ 



sf3 2 -c + iO -s(3 2 - d + iO 



-1, (8) 



and a similar integral over the variable ct\. After all operations we can write the matrix 
element in the form 

{2) _ . (Wtc^Z.Z^N.N, f d 2 fc!d 2 fc 2 u(q-)Rv(q + ) 

(2) ^ (2!) 2 J vr 2 klkl{ qi -k,y{q 2 -k 2 r [) 

with 

B. The classification of diagrams 

It is convenient to classify FD by order of exchanged photons absorbed by the lepton 
world line (Fig. |2J). We mark them as Rijki, R = X] Rijki with different integers i, j, k, I from 
one to four, counting from a negative lepton emission point. 

a) Consider first the set of 4 FD (Fig. 0k) named -R1234, -R2134, -R1243, -R2143 in which 
the interactions with two nuclei are ordered consecutively against the lepton line direction. 
The sum of relevant contributions provides the convergence of {3i,a 2 integrations. After a 
standard calculation one obtains for this set 



-Rl234 + -^2134 + -^1243 + -^2143 — "5 2 — , a 1 v? — — —B — , 



q. ± (q. -q 1 ) ± 

B ~ r<r • >i<i <ir (10) 

The last equality in (fTU|) is the result of Dirac equation for massless particles 

u(q-)P-pip 2 = -u(q-)q-±p 2 . (11) 
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A similar result as in (fTUj) is achieved for the set of crossing diagrams (Fig. HJd) relevant to 
-R3412, -R3421, -R4312, -R4321 terms in the amplitude with only the replacement B — > B where 
for B stands 

5_ (~g+ + gi)±g+± / 19 x 
/3_<£ + M<?i-<?+) 2 ' 1 J 

b) Let us now consider the set of the diagrams -R1342, -R1432, -R23415 -R243i(Fig. Eh) and 
-R3124, -R3214, -R4123, -R4213 (Fig. Hi), where exchanges with the ion B(A) are attached to the 
lepton line between the interaction with the ion A(B). 

For defmiteness, consider the sum -R1342 + -^1432- Using the relevant denominators of the 
lepton line one obtains the following integrals over /3±, a 2 - 

d/3i 1 

2vri sa.(f3_ - fa) - (q_ - k{) 2 + iO 
1 



x 



x 



+ 



-sa + {(3_ - Pi) - {-q+ + qi - k,) 2 + iO 



da 2 

27Tl" 



s(J3- - A) 



s(/L - A)(a_ - a 2 ) - (q_ - k Y ) 2 + iO 
s(P--Pi) 



(13) 



s{(3_ - (3i)(-a+ + a 2 ) - (-q + + q 1 - k ± ) 2 + iO. 

The second integral after closing the integration contour in the lower half plane gives the 
function sgn(/3_ — /3%), thus (fTB^) becomes 

dfa sgn(/3i-/3_) 



(14) 



2vri sa.{/3_ - ft) - (q_ - k^ 2 + iO 
1 

X (sa+(/3- - K) ~ (-<?+ + q~i ~ fci) 2 + iO) ' 
Using the relation 

no 

dx sgn(x) 1 , ad 

2ni (— ax — b + iO)(cx — d + iO) iri (ad + be) be 1 
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we obtained the following result: 



-R1342 + -R1432 + -R2341 + -R2431 



Pi f (q--ki)±(-q+ + qi-ki)± a+(q_-ki) 



ins \a + (q_ - fa) 2 + + q 1 - k ± ) 2 q + + q 1 - ki) 2 

| (g- ~ gi + fa)±(-q+ + fa)± ^ a + (q_-g 1 + fa) 2 \ 
a + (q_-q 1 + k 1 ) 2 + a-(-q + + k 1 ) 2 U a.(-q + + fa) 2 / ' 



-R3124 + -R3214 + -R4123 + -R4213 



fa ( (q--fa)±(-q+ + q2-fa)± , P-(-q+ + <? 2 - k ^ 2 



ms \(3 + (q_ - k 2 ) 2 + /?_(-«+ + <? 2 - fc 2 ) 2 ln ' (3 + (q_ - k 2 ) 2 



{q- - q_i + k 2 )±{-q+ + k 2 )± . (3-(-q + + k 



\2 



0+(q- -Q2 + fa) 2 + 0-(-q + + k 2 ) 2 11 - q 2 + k 2 ) 2 ) ' 1 ' 

It is necessary to point out that the obtained expressions are pure imaginary and 

consequently their interference with the Born term in the cross section is zero. 

c) Consider the case of interactions with different nuclei alternating along the lepton line, 
for instance, the amplitude -R1324 (Fig. 0^). After a bit algebra one obtains for the relevant 
numerator 

A^i324 = sf>ifi 2 {q„ - h)±(q- -fa- k 2 )±(q- -q\ - k 2 )±, (17) 

which is very different from the numerators of the Born like amplitudes. Specifically, it is 
the term of higher order in the running transverse momenta ki. 
The relevant denominators read 

{1} = (g_ - fa) 2 + i0 = s(/L - ft)a_ - (q_ - fa) 2 + iO, (18) 
{2} = (q_ -fa- fa) 2 + iO = s(/3_ - ft)(a_ - a 2 ) - (q_ -fa- k 2 ) 2 + iO, 
{3} = (-<?+ + g 2 - fa) 2 + iO = s(-/3+)(a_ - a 2 ) - (-q + + q 2 - k 2 ) 2 + iO. 

The nonvanishing contribution only emerges if the poles are located in different a 2 half- 
planes, which takes place only if p 1 < j3_ (/3± > 0). Taking the residue at the pole {2} we 
find 

sda 2 1 0(0- -ft) (ig) 



2m {2}{3} (P 1 -[3_)(-q + + q 2 -k 2 ) 2 -l3 + (q_-fa-k 2 ) 2 - 
Further integration over j3\ can be done using the relation 

run 

dx 9{x) 1 ^i+iiB^y (20) 



27ri (ax - b + i0)(cx + d + iO) 2(ad + bc)\ n be 
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wit the result 

/3-iV"i324 A , i , ad\ 



^1324 = -^^ 1 + -In- , 
2sDi324 V n be J 

Dn24 = P-{q- - k 1 f{-q + + q 2 - k 2 f + (3 + q 2 _ (q_ - k x - k 2 ) 2 = ad + be. (21) 

The highly nonlinear denominator (|2*TJl makes the contribution from the considered case 
dramatically different from the Born like amplitude. Technically, the nonlinearity is not 
surprising because of the related nonlinearity of the numerator. The principal difference 
from the Born like amplitude is that with the alternating ordering of interactions we have 
the situation in which the p + component of the lightcone momentum is conserved in the 
scattering on one ion but is not conserved in the scattering on the second ion. Depending 
on the ordering of interaction vertices and the order of integrations one will encounter 
the mismatch of conservation and nonconservation of the p- component of the lightcone 
momentum. 

Similar results can be obtained for other contributions of these types, 
d) The final result reads 



<U = W (™™ ^ — — fcjtfo, -»,)»(«, -*,)' • 



,(2) [a>nb n }± 



To convince the gauge invariance fulfilment we put the explicit form for the real part of 
the amplitude 

Ro d(2) _ [g-(g- -gi)U , [(-9+ + Qi)q+}± 



[g-(g- - fciKg- - fci - fcgKg- - gi - 

2[/3_(q_ - + q 2 - fc 2 )2 + (3 + ql(q_ - - k 2 f] 
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n 






K 






1 


R(12)(34) 




Q- 


q- - qi 






2 




Qi 


— q + 


q+ 






3 


#1324 




Q- 


q- - h 


q- — k\ — k 2 


Q-—Q1 — k 2 


4 


#1423 




Q- 


q- - h 


q- - q2 + k 2 - ki 


-q + + k 2 


5 


#2314 




<7- 


Q--Q1 + h 


q- - qi + h - k 2 


-q+ + 92 - k 2 


6 


#2413 




q- 


q- - qi + h 


—q + + k x + k 2 


-q+ + k 2 


7 


#4231 


q- - 


<72 + h 


—q+ + ki + k 2 


-q + + k\ 


q+ 


8 


#3241 


Q— 


- k 2 


q- - qi + h - k 2 


-q + + k\ 


q+ 


9 


#4132 


q- - 


<72 + h 


q- — <?2 + k 2 — k\ 


—q + +q\ — k\ 


q+ 


10 


#3142 


Q- 


-k 2 


q- — k\ — k 2 


-q+ + qi- h 


q+ 


11 


#3(12)4 


Q- 


-k 2 


-q + +q 2 -k 2 






12 


#4(12)3 


q- - 


<72 + k 2 


-q+ + k 2 







Tab. I: The coefficients for formula Q23|) . The brackets denote index permutation, e. g., (12) 
12 ■ 21. 

[g-(g- - fci)(g- -q2 + k 2 - fci)(-g+ + k 2 )}± 

2[P4q_ - fci) 2 (- 9+ + k 2 f + (3 + ql(q_ ~q 2 + k 2 - k^] 

[q-(q- -q\ + fciKg- - gi + h - k 2 )(-q + + q 2 - k 2 )} ± 

~ 2[(34q_ -q, + k x f{-q + + q 2 - k 2 f + {3 + q 2 _(q_ -q 1 + k 1 - k 2 f] 

[q-{q- -qi + k 1 )(-q + + h + k 2 )(-q + + k 2 )]± 

2[P-(q- -Qi + ki) 2 (-q + + k 2 f + (3 + ql{-q + + k x + k 2 f] 

Kg- - ?2 + k 2 )(-g+ + h + fc 2 )(-g+ + fci)g+]± 

2[/9_^(-q + + fci + k 2 f + /?+(-<?+ + fci) 2 ^ - q 2 + fc 2 ) 2 ] 

[(g- ~ fcgKg- - gi + h - k 2 )(-q + + k 1 )q + ]± 

~ 2[f3_q\{q_ -q x + k x - k 2 f + (5 + {-q + + k,)^ - k 2 f] 

Kg- - 92 + k 2 )(g- -q 2 + k 2 - k 1 )(-g+ +_gi_- fci)g+]± 

2[/3_q 2 h (<7_ ~q 2 + k 2 - k,f + - <? 2 + fc 2 ) 2 (-q + + 9l - A^) 2 ] 

Kg- - k 2 )(q- -h - k 2 )(-q + +_gi_- fcpg+jx 

2[/3-9 2 h (9_ - fci - k 2 f + (3 + (-q + + k,)^ - k 2 f] " 

Then one can verify that the following condition is satisfied: 

Rei?[ 2 j=0 if fei = or k 2 = or k ± = q x or k 2 = q 2 . (24) 
This fact is correct also for the whole amplitude (J23j). As one can see, this property is crucial 
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for the infrared convergence in integrations over d 2 ki . 

Under the loop integration one can make the shift of the integration variable ki — > q i — ki. 
Then expression (J2H|) for Re R,J can be simplified to 

R d(2) g-±(g- ~Ql)± , (~g+ + 9l)±9+± 

™«P) ~ M _ + p _ {q _ _ ^2 + p _ q 2 + + ^ _ g+) 2 

_ 2 - fciKg- - fci - fegKg- - gi - h)}± 

(3-{q_ - fci) 2 (-q + + q 2 - fc 2 ) 2 + /3 + g 2 (q_ - kt - fc 2 ) 2 

[(-g+ + gi + fc2)(-g+ + fci + fc2)(-g+ + h)q+}± 
P_ql(-q + + k 1 + k 2 f + (3 + (-q + + k i y(q_-q 2 + k 2 y 1 J 

Despite the gauge invariance property is not seen clearly her,e as in the previous case, the 
final results after integration over ki coincide. 

III. THE WIDE ANGLE LIMIT OF THE Af ( S AMPLITUDE 

Let us consider the behavior of this expression in the case when the transverse component 
of lepton momenta is large compared to the momenta transferred to the ions 

<?- ~ -q+ = q, \q\ > l<?i, 2 |- (26) 

In this case, the main contribution to the matrix element gives the region 

Wi\ < \H < \q\. (27) 



The amplitude M^j reads 



,(i) 



= _ i3 (sM^ ^ v(q+): (2g) 

H1H2 b 

1 * 9- - 9i * , * 9i - 9+ . / n „\ - 



For wide angle kinematics one has 



1 ^ 1 

- Ri (X) = -irn^^q^mi + 2/3-q. qi ] + q 2 [b iqi q 2 + 2(3 +qi .q 2 \], (29) 



s y ' s 



bl(q 



with b\ = P_ + P + , q = q ~ — q + and q 1 2 are the transferred to ions momenta. 
For matrix element M^) we have (in agreement with the result obtained in paper 

(1) _ 2W Zl zl Nl N 2 r <Pk 2 i )R wPi v(a ) m 



.is]; 
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with 



,(i) _ D , s (q--k 2 )±(q--qi-k 2 )± (q+ - h - q 1 ) ± (q+ - k 2 ) 



R{2) B + B /?-(<?- ~ <h - k 2 y + (3 + (q_ - k 2 y (3 + (q_ - q 2 + k 2 f + ^(q + - k 2 f 

(31) 

In the considered limit this expression has the form 



R {1) 



l2 > b lQ 2 

_ 2g,.fc 2 
g 2 



(2/3_g_.g 1 + g_g 1 )( 



4( g_.fc 2 ) 2 k\ 
{q2)2 q2 



(k 2 qi + W-k 2 .q x )\ + ((3. -> |fe 2 | > |g 2 |. (32) 



This expression turns to zero after the angular averaging. It can be shown that the quantity 
M$ as well turns to zero in the limit of wide angles pair production and is proportional to 
\q 2 1/| q\ <C 1, which is in agreement with 

For the considered above amplitude M$ (l2^j) the quantity Rfl plays a role of cut-off in 



'(2) yany vjuc*j.j.uj.ujr ^t-/ 2 ) 
the region |fej| > |qr|. From very general arguments it can be cast in the form 

Rp R (2) „ [Kiqi-hYk^-hrU 

rve n,^2) ~ (<J 2 ) 2 ^LivaPi \°°) 

with some dimensionless tensor matrix independent of fej, q v Expanding the expres- 

sion (J25J) one gets 

f d^k 2 Reggj ^ I 4(/3 + -/3_) q 2 <?L. 

J vr 2 k\kl{ qi - fe 1 ) 2 (q 2 - fc 2 ) 2 ~ (<Z 2 ) 2 (P- + P + ) 2 q\ q\ ' 1 J 

where J is the unit matrix and q max is the upper integration limit q max — -R is the 

nucleus radius. Such enhancement is absent if the number of exchanged photons from every 
ion exceeds two (Fig. EJ). Really, the amplitudes m9\, M$ , n > 2 contain only the first 



power of large logarithm, whereas m, n > 2 do not contain such a factor at all, because 

the corresponding loop momenta integrals are convergent in both infrared and ultraviolet 
regions and one can safely put |<7x(2)l = over loop integrations. 

Thus, the differential cross section for the considered kinematics is determined by the 
interference term (M^j)* Mm! which has the form (for comparison we present also the Born 
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term) 



d„ _ 1 6 ( w)M^ + (i-xn d2?id2?A 



dfoidx vr 4 q\q\{q 2 ) 2 h 

da mt _ IGjZ^a 2 ) 3 (1 - 2s) ^ gg^. ^ g 2 
dfoidx q\q\q\q 2 - &i q\ q 2 



l n ^^l n ^^gd 2 g 1 d 2 g 2 d 2 g_ ) (36) 
"2 



- x)q 2 _ + x(g_ - q^ 2 xq\ + (1 - x)^ - <? + ) 2 ' 

/3_ , 1 4m 2 x(l — x) 
x = — , e < sc, Oi < 1 — e, e = ^ 



We note that expression (J3HJ) is symmetric under simultaneous substitutions g + <->• g_ and 
/3 + <-> /3_ due to the C-even character of the interference. 

Finally, from very straightforward generalization of (}3*3"j) it can be shown that the set of 
amplitudes with an odd number of exchanges with one or both nuclei is suppressed in the 
limit of wide angle production 

*C, r o(^), ^sr'-Kit/)' ~°(n t ^f d )' (37) 

IV. MULTIPHOTON EXCHANGE 

Let us generalize the above picture for the case of multiple photon exchanges (m, n > 2). 
Using the relation 



j 1 f d 2 fcx ■ ■ ■ cgfcg-i _ nln"- 1 (^-) 

» Trn-iy (fe 2 + A 2 )...(fe 2 _ 1 + A 2 )((q-fe 1 fc n _ 1 )2 + A 2)- q 2 ' 1 ) 

and taking into account the combinatorial factor -^ coming from the symmetric integration 
over «j, one has to replace any single photon exchange by an infinite set of photons, 
multiplying the amplitude by the factors of type exp{k/?j(qr 2 )} with the phase <fii(q 2 ) = 
±aZi In The scattering of electron and positron differs only by sign of the phase (positive 
for electrons) [^]. This replacement is depicted in Fig. |H1 where the double photon line 
corresponds to the infinite set of photons. 

Using the same technique as in one can see that the amplitude relevant to Fig. [T^l 
and Fig. [TJd can be cast in the form 

= Be -i[<fi(Qi) ~ <f2(ql)) + se^tei) ~ vMl)}. (39) 
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The interaction of the electron and the positron with Coulomb field differs only by signs. 
Though this expression is infrared unstable in the case Z\ 7^ Z 2 the regularization parameter 
A enters it in a standard way. 



Let us now consider the class of diagrams depicted on Fig. [Zb- In subsection III B| we 
obtained the expressions pfia[ I16bjl for the case m = n = 2 such that Rei?i( 34 ) 2 = 0. It can 
be shown that the terms of higher order with any even number of photons from same nuclei 
attached to the lepton world line between two photons from other nuclei do not contribute 
to the amplitude of the process under consideration. It is the consequence of the relation 
(sgn(a)) 2fe+1 = sgn(a). 

The general structure of the amplitude corresponding to Fig. [2b can be constructed using 
the lowest order truncated amplitude (without single photon propagators) R^ 1 ' 



(2) 



COS^gf) ( 1} iMfc 2 } _ (( _ fe)2)] 
^(2) ^2 rL {2)* 3 



= 1 (q--q 2 + k) ± (-q + + k) ± P+(q--q 2 + k) 2 

L ( 2 ) i7rf3_( q+ -k)i + f3 + (q_-q 2 + ky ° P-{q + -kf 



The further generalization is obvious. For instance, we cite the expression corresponding 
to the diagram depicted on Fig. 01 

R$ = cos^fc^e-Wfai " fc i) 2 ) cos^ 2 (^)e i ^((92 " k 2 f) R ^ (41) 

From the above consideration we conclude that the general structure of the matrix element 
, corresponding to m photon exchanges from one ion and n exchanges from other one, 
schematically reads 



r (m) • ~ ~ v 



2 



M"> = isNMZ^) m (Z 2 OL 



( ' 16n!m! 

x j 1 1 ^^ )t (42) 



where m and n obey the condition \m — n\ < 1. At this stage, we omitted phase factors in 

,(m) 
L (n) 



the structure (for clearly understanding the problem), so it can be written in the form 



5(m) _ 5(1) 5(1) 5(2) 5(2) 5(2) 5(3) p(3)i? p(3)L ,,„v 
K (n) ~ K (l) + ^(2) + ^(1) + K {2) + ^(3) + K {2) + ^(3) + ^(3) " ■ ■ l 4c V 



15 



R 



(2) 
(1) 



1 {q--k)±(-q+ + qi-k)_ 



R {2) 



i7r a_ 



R i3) 



q+ + Qi ~k) 2 + a + (q„ - k)'- 



In 



a + {q_ 



k) 



OC- 



{-Q+ + Qi - k) 



2 ' 



0. 



0)R _ 1 
( 3 ) ~ Ci + c 2 



2 c l 



2 c 2 



c 2 



7T 



Cl 
C 2 



/?_(<?_ - k 1 f{q_ - fci - fc 2 - Ki) 2 (-g + + q 2 - Ki - K 2 ) s 



(3) 
('1) 



= P+Q-(Q- - ki - K 1 ) 2 (q_ - fci - k 2 - «i 

— — 3( 2 + - In — 

di + d 2 L 2 a 2 J 



«2 



^1 = P+(Q- - K l) 2 (^- - Kl — K 2 - ^l) 2 (q- - Ki - K 2 ~ K 3 - k X - fc 2 ) 2 , 

rf 2 = - Ki- k 1 ) 2 (q_ - m - K 2 - fci - fc 2 ) 2 (-qr + + q 2 - kj - k 2 - k 3 ) 2 . (44) 

Here i?| 2 j is only the second term in the right-hand side in (f2"H|) and the index R(L) denotes 
two possible configurations of photons for R^ R (Fig. Ek) and R^ L (Fig. [Tf). 

In such a way, the general algorithm for construction of an arbitrary term is transparent. 
Unfortunately, we cannot obtain the compact expression for the whole amplitude. The 
reason is the increasing nonlinearity of the propagators with the order of interaction. The 
behavior of the above denominators is very different from the Born-like case, where the 
simplicity of propagators allows one to obtain eikonal-like expressions. 



V. CONCLUSIONS 



The wide angle lepton pair production in peripheral interactions of ultrarelativistic heavy 
ions is an archetype reaction for hard processes in central hadronic hard collisions of heavy 
nuclei. In the electromagnetic case, the expansion parameter Z 12 a ~ 1 makes the multiple 
photon collisions, m^ + n^ — ► 1 + U potentially important ones, likewise the effect of multiple 
gluon collisions in central collisions is enhanced by a large number of nucleons at the same 
impact parameter. The crucial issue is whether such multiple photon collisions can be 
described by the Born cross section in terms of the collective photon fields of colliding nuclei 
or not. We obtained the expression for the amplitude 27 + 27 — * 1 + 1~ and show that its 
contribution is dominant in a wide angle limit. Our principal finding is that the amplitude 
is manifestly of non-Born nature, which is suggestive of the complete failure of linear k±- 
factorization even in the Abelian case. 
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We have shown that the terms in perturbation series of the amplitude for the process of 
lepton pair production in the Coulomb fields of two relativistic nuclei relevant to the closed 
two photon loops are logarithmically enhanced in this case, while in higher order terms 
such enhancement is absent. We presented the algorithm which allows one to construct the 
full amplitude in all orders. The obtained results can be useful in application to the QCD 
process of production of high kj_ jets, the issue which will be investigated elsewhere. 
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Fig. 1: Typical Feynman diagram for the amplitude M 
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(2) 




1 n 
+ all permutations 



Fig. 2: The notation for the permutations of n virtual photons emitted by the heavy ion. 
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Fig. 3: The set of basic Feynman diagrams for the amplitude M 
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Fig. 7: The Feynman diagrams for the amplitudes with many photon exchanges. The double 
photon line represents any number of exchanged photons, the double zigzag line represents only 
the odd number of exchanged photons. 



